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ABSTRACT 

We present measurements of the dark matter bispectrum in N-body simulations with 

non-Gaussian initial conditions of the local kind for a large variety of triangular con- 
figurations and compare them with predictions from Eulerian Perturbation Theory up 
to one-loop corrections. We find that the effects of primordial non-Gaussianity at large 
scales, when compared to Perturbation Theory, are well described by the initial com- 
ponent of the matter bispectrum, linearly extrapolated at the redshift of interest. In 
addition, wc find that, for Jjsil = 100, the nonlinear corrections due to non-Gaussian 
initial conditions are of the order of ^ 3-4% for generic triangles up to ~ 20% for 
squeezed configurations, at any redshift. We show that the predictions of Perturbation 
Theory at tree-level fail to describe the simulation results at redshift ^ = at scales 
corresponding to fc ~ 0.02-0.08 /iMpc~^, depending on the triangle, while one-loop 
corrections can significantly extend their validity to smaller scales. At higher redshift, 
one-loop Perturbation Theory provides indeed quite accurate predictions, particularly 
with respect to the relative correction due to primordial non-Gaussianity. 
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1 INTRODUCTION 

In recent years, a significant research activity has been 
devoted to the effects of a possible small departure from 
Gaussianity in the primordial cosmological perturbations. 
While current constraints on primordial non-Gaussianity 
from measurements of the cosmic microwave background 
(CMB) and the large-scale structure are still consistent 
with the Gaussian hypothesis (Komatsu et al. 2010; Slosar 
ct al. 2008), a possible detection in forthcoming experiments 
would constitute a major discovery, providing crucial infor- 
mation on the early Universe and on the high-energy physics 
of inflation (sec, for instance, Komatsu et al. 2009). 

The effect of primordial non-Gaussianity on the large- 
scale structure has been assumed, for a long time, to be 
limited to an additional, primordial component to the mat- 
ter skewness and bispectrum induced by gravitational in- 
stability and to a correction to the abundance of massive 
cluster (see Liguori et al. 2010; Desjacques & Scljak 2010b, 
for recent reviews). Numerical and analytical studies have 
indeed shown that a matter density probability distribution 
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initially skewed toward positive values produces more over- 
dense regions and, consequently, collapsed objects while a 
negatively skewed distribution produces larger voids (see 
Grossi et al. 2008; Pillepich et al. 2010; Kamionkowski et al. 
2009; Lam & Sheth 2009; Lam et al. 2009; Maggiore & Ri- 
otto 2009, for recent work). Moreover, a nonvanisliing skew- 
ness in the initial conditions corresponds to a primordial 
component to the matter bispectrum, i.e. the three-point 
function in Fourier space. For the local non-Gaussian model 
considered here, the primordial matter bispectrum exhibits 
a scale, redshift and triangle shape dependence distinct from 
that of tlic component sourced by the nonlinear growth of 
structures. This enables us in principle to disentangle the 
two contributions. In the specific case of equilateral triangu- 
lar configurations, the primordial contribution to the matter 
bispectrum scales as ~ kr'^ relative to the gravity-induced 
term, leading to large, potentially observable corrections at 
low wavenumbers. Measurements of the galaxy bispectrum 
in future large-volume redshift surveys (such as Euclid or 
HETDEX) should be able to provide constraints on the lo- 
cal non-Gaussian model competitive with those from CMB 
observations (Scoccimarro et al. 2004; Sefusatti & Komatsu 
2007; Sefusatti et al. 2009). 

In addition to these effects, Dalai et al. (2008) have re- 
cently discovered a large correction to the galaxy bias in 
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numerical simulations of local primordial non-Gaussianity. 
Further numerical and theoretical work has confirmed this 
result (Matarrese & Verde 2008; Slosar et al. 2008; Afshordi 
& ToUey 2008; McDonald 2008; Grossi et al. 2008; Taruya 
et al. 2008; Pillepich et al. 2010; Desjacquos et al. 2009; 
Giannantonio & Porciani 2009). The constraints obtained 
from power spectrum measurement of highly biased objects 
in current datsrsets are already comparable to the CMB 
results (Slosar et al. 2008; Desjacques & Seljak 2010a), and 
the prospects for detecting local primordial non-Gaussianity 
with galaxy clustering look exciting (Dalai et al. 2008; Car- 
bone et al. 2008; Seljak 2009; Slosar 2009; Verde & Matarrese 
2009; Desjacques & Seljak 2010a). At this point, analyses of 
the galaxy bispectrum preceeding the work of Dalai et al. 
(2008) must be updated to account for the non-Gaussian 
correction to the galaxy bias. In fact, a rigorous joint analy- 
sis of the galaxy power spectrum and bispectrum in presence 
of local non-Gaussianity is in order. First steps in this di- 
rection have been taken by Joong & Komatsu (2009) and 
Sefusatti (2009) with a preliminary comparison with simu- 
lations in Nishimichi et al. (2009). 

In this perspective, we will consider the measurement 
of several triangular configurations of the matter bispectrum 
on mildly nonlinear scales, with both Gaussian and non- 
Gaussian initial conditions of the local type. Although the 
matter bispectrum is not directly observable with tracers 
of the large-scale structure, it is instructive to assess the 
extent to which perturbation theory describes the shape de- 
pendence of the matter three-point function in the presence 
of non-Gaussianity of the local type. This analysis will be 
useful when considering the complication brought by bias- 
ing, which will be addressed in a forthcoming publication. 
Measurements of the matter power spectrum with local non- 
Gaussianity can be found in Pillepich et al. (2010); Des- 
jacques et al. (2009), where the small corrections at mildly 
nonlinear scales predicted in the framework of perturbation 
theory by Taruya et al. (2008) are observed. In the case of 
the matter bispectrum, measurements in simulations with 
Gaussian initial conditions are shown in Scoccimarro et al. 
(1998); Hou et al. (2005); Pan et al. (2007); Smith et al. 
(2008); Guo & Jing (2009), with Smith et al. (2008) con- 
sidering, in addition, rcdsliift space predictions in the con- 
text of the halo model. By contrast, the only measurement 
so far of the matter (and halo) bispectrum in simulations 
with local non-Gaussian initial conditions can be found in 
Nishimichi et al. (2009), where a relatively small subset of 
isosceles triangular configurations is considered. 

We will compare our measurements with predictions 
of the matter bispectrum at the one-loop approximation in 
Eulerian perturbation theory. A comparison of one-loop re- 
sults with the bispectrum extracted from simulations with 
Gaussian initial conditions is shown in Scoccimarro et al. 
(1998), whereas a comparison of the effect of primordial 
non-Gaussianity with the tree-level prediction of pertur- 
bation theory is performed in Nishimichi et al. (2009) for 
"squeezed" isosceles configurations &t z — with k < 
O.l/iMpc"^ only. Here, we will extend the analysis to in- 
clude several triangular configurations covering the range of 
scales 0.002 < A; < 0.3/iMpc"^ and redshifts z = 0, 1 and 
2. This will allow us to broadly test the accuracy of one- 
loop perturbation theory in the mildly nonlinear regime. We 
will also discuss the validity of two phenomenological pre- 



scriptions for the nonlinear bispectrum with Gaussian ini- 
tial conditions, namely the fitting function of Scoccimarro 
& Gouchman (2001) and the formula of Pan et al. (2007) 
based on a scaling transformation. 

This paper is organized as follows. In section 2 we sum- 
marize previous results on the predictions of the matter 
power spectrum and bispectrum in cosmological perturba- 
tion theory for both Gaussian and local non-Gaussian initial 
perturbations. In section 3 we describe the N-body simula- 
tions and the bispectrum estimator employed in our analysis 
whereas, in section 4, we present our measurements of the 
matter bispectrum and compare them to one-loop predic- 
tions in perturbation theory. Finally, we conclude in sec- 
tion 5. 



2 THEORY 

In this section, we summarize previous results on the nonlin- 
ear evolution of the matter correlators as described specifi- 
cally by Eulerian Perturbation Theory (PT). The quantity 
of interest, the matter overdensity S, is obtained as a per- 
turbative solution to the continuity and Euler equations, 
and Poisson equation relating matter perturbations and the 
gravitational potential. These equations fully determine the 
evolution of the matter density and velocity fields, once the 
initial conditions are given in terms of the primordial cor- 
relators. Other approaches such as Lagrangian Perturbation 
Theory, for instance, have also been studied in the literature. 
We refer the reader to Scoccimarro (2000) for a study of the 
matter bispectrum in Lagrangian Perturbation Theory with 
Gaussian initial conditions and to Bernardeau et al. (2002) 
for a complete review of cosmological perturbation theory. 

2.1 Initial conditions 

Our N-body simulations of the matter density evolution as- 
sume local non- Gaussian initial conditions. This model of 
primordial non-Gaussianity is defined by the local expression 
in position space for the Bardeen's curvature perturbations 
$ (Salopek & Bond 1990, 1991; Gangui et al. 1994; Verde 
et al. 2000; Komatsu & Spergel 2001) 

$(x) = <^(x) + fr,L [0'(x) - (<^2(x))] , (1) 

where the second term on the r.h.s. represents a non- 
Gaussian correction to the Gaussian random field </>(x). In 
this expression, we assume that ^(x) is the curvature field 
during early matter domination, and not the linearly extrap- 
olated value at present time. Despite its relatively simple 
form, the parameterization of primordial non-Gaussianity 
provided by equation (1) well describes inflationary models 
in which the non-Gaussianity is produced by local mecha- 
nisms on superhorizon scales (see Bartolo et al. 2004; Liguori 
et al. 2010; Chen 2010; Byrnes & Choi 2010, and references 
therein) . 

The definition of equation (1) corresponds to a very 
specific functional form of the bispectrum and trispectrum of 
the initial curvature perturbations. One finds the following 
leading contribution to the curvature bispectrum, 

Ri>(fci, fc2, hi) = 2/ivLP#(fci)P$(fc2) + 2 perm. , (2) 
with the curvature power spectrum P#(A:) defined in terms 



of the Gaussian component alone as {(f){'ki)(j)(k2)) = 
(5^^(ki2)P*(fci), where we introduce the notation ktj = 
kj + kj. The magnitude of the curvature bispectrum is maxi- 
mized for "squeezed" triangular configuration, i.e. when one 

side of the triangle is much smaller than the other two, say 
fci ^ ^2 — fea- The curvature trispcctrum is given by, 

T4.(ki,k2,k3,k4) = 4/^LP<i.(fci)P#(fc2) X 

[P^ikia) + P*(fci4)] + 5 perm. (3) 

The linear matter overdensity in Fourier space (5k is 
related to the curvature perturbations $k via the Poisson 

equation. 



<5k(a) = M(fc,«) $k, 
where we introduced the function 

2 k'^T{k)D{z) 



M{k,z) 



3 ^0 



(4) 



(5) 



with T{k) being the matter transfer function, computed with 
the CAMB code (Lewis et al. 2000), and D{z) the growth fac- 
tor in units of 1 + z. The initial matter correlators arc related 
to the correlators of the curvature perturbations through 

(5ki • • • ^k„) = M{ki,z) ■ ■ ■ M{k„, z){$ki ■ ■ ■ $k„) , (6) 

so that the linear power spectrum is given by 

Po(k) = M''{k,z)Pi,{k), (7) 

while the initial bispectrum and trispectrum are given re- 
spectively by 

Bo{ki,k2, ks) = M{ki)M{k2)M{k3)Bi{ku fe, ks) , (8) 
and 

ro(ki,k2,k3,k4) = M(ki)M{k2)M{kz)M{ki) X 

r<l.(kl,k2,k3,k4). (9) 

As we will see shortly, nonlinear corrections to the matter 
bispectrum will depend on both the initial bispectrum Bo 
and trispectrum Tq. 

2.2 Perturbation theory 

In PT, the solution for the nonlinear matter density contrast 
Jk in Fourier space is given in terms of corrections to the 
linear solution 5^^^ (Fry 1984), so that 



where each nonlinear correction given by 



(5*"' = j (fqi . . . d'^q„Fn{q_i, . . . , q„ 



(10) 



(11) 



with F„(qi, . . . ,q„) the symmetrized kernel of the n-order 
solution. Equation (10) allows one to derive the evolved mat- 
ter correlators once the initial correlators, i.e. the correlators 
of the linear are known. For Gaussian initial conditions, 
only the linear power spectrum Pq must be specified. In gen- 
eral however, higher-order correlators need to be taken into 
account. 

In analogy with quantum field theory, pcrturbativc so- 
lutions for the matter correlators can be denoted as tree-level 
or one-loop, two-loop, etc., according to the number of in- 
ternal integrations present in their expressions. However, it 



should be noted that, while in the case of Gaussian initial 
conditions the number of loops of the perturbative correc- 
tion correspond univocally to a specific perturbative order, 
this is, as we will see below, no longer true for non-Gaussian 
initial conditions. 

For completeness, we summarize here the explicit ex- 
pressions of the one-loop PT expansion for both the matter 
power spectrum and bispectrum with generic non-Gaussian 
initial conditions. In the case of the matter power spectrum, 
we have up to fourth order in PT (see Bernardeau et al. 
2002, and reference therein) 



P{k) = Pll(fc) + Pl2(fe) + Pi2 + Pis + 

two— loop terms + 0{5q), 



(12) 



where, Pn = Po is the linear matter power spectrum, while 



Pi2 = 2 /d3gP2(q,k-q) Bo(fc,g,|k-q|), 



(13) 

P2'2 = 2 Jd'qF^iq, k - q) Po{q) Po(|k - q|), (14) 

P/a = 6Po(fc)y"dV3(k,q,-q) Po(9). (15) 

We can see that the only additional contribution due to pri- 
mordial non-Gaussianity is Pi2(fc) which depends on the ini- 
tial bispectrum Bq (neglecting two-loop contributions at the 
fourth order in PT that depend on the initial trispectrum). 
The amplitude of this correction for local non-Gaussian ini- 
tial conditions was studied in Taruya et al. (2008), who con- 
sidered also initial conditions of the equilateral kind. They 
found that the effect of a primordial non-Gaussian com- 
ponent within the bounds from CMB observations is typ- 
ically below 1% at mildly nonlinear scales, at the limit of 
detectability in future large-scale structure observations. 

One-loop corrections to the matter bispectrum for 
Gaussian initial conditions have been studied in Scoccimarro 
(1997); Scoccimarro et al. (1998) while the extension to 
generic non-Gaussian initial conditions is explored in Se- 
fusatti (2009). For the bispectrum up to sixth-order in PT 
and excluding two-loop corrections, we have the following 
expression 

B = -Bill -|- Bii2 + ^112 + 

Pi22 + Pi22 + Piia + Pns + 

P222 + P123 + Pila + Pii4 + 2-loop terms, (16) 
where Sm s Bo is the initial bispectrum and 

B(i2 = 2 P2(ki,k2) Po(A;i) Po(fc2) + 2 perm., (17) 
is the other tree-level contribution, while the 1-loop correc- 



tions are given by 



■t>112 



- P 



q F2(q,k3-q) To (ki , k2 , q, ks - q) 



(18) 



2 Po(A;i) 



F2 (ki , ks) i\ F2 (q, kg - q) x 



B 



Biis 



-B222 

Bi OS 



-t>l23 



Bo{k3, q, |ks - q|) + (fcs -H- fe)] + 2 perm. 
F2(ki, k2) [Po(A;i) Pi2(fe2) + Po{k2) Pi2{ki)] + 
2 perm., (19) 

4 y"rf^g F2(q,k2-q) F2(ki +q, ka-q) x 

Bo(fci,g, |ki+ql) Po(|k2-q|) 

+2 perm., (20) 

3Bo(fci,fc2,fc3) y^d^ P3(k3,q, -q)Po(g) + 

2 perm., (21) 



3Po (fci )y d g Ps (ki , q, ka - q)Bo (fe , g, I k2 - q| ) + 
(fci o fe) + 2 perm., (22) 
8yrf''gP2(-q,q+ki)P2(-q-ki,q-k2) x 
P2(k2-q,q)Po(g)Po(|ki+q|)Po(|k2-q|), (23) 
6 Po(fci)yd^g P3(ki,k2-q,q) P2(k2-q,q) x 
Po(|k2-q|) Po{q) + 5 perm., (24) 

6 Po(fci) P„(/C2) P2(ki,k2) X 

J df'q Pi(ki, q, -q) Po((j) + 5 perm. 

P2(kl,k2) [Po(fcl) Pl3(fc2) + Po(fc2) Pl3(fcl)] + 

2 perm., (25) 
12 Po(fci) Po(A;2)yd^gP4(q,-q,-ki,-k2) Po(9) + 
2 perm.. (26) 



Specifically, the one-loop contributions present because of 

non-Gaussian initial conditions are B112, which depends on 
the initial trispectrum To, and all the fifth-order terms -6122, 
^122 1 ^113 and Sffs, which depend on the initial bispectrum 
So. 

The remaining terms, corresponding to Gaussian initial 

conditions, were recently studied in the context of rcsum- 
mation techniques of the PT series and can be regarded as 
perturbative expansions of "resummed" kernels (Bcrnardeau 
et al. 2008). For instance P(|s corresponds to the next-to- 
leading term in the resummation of the nonlinear propa- 
gator in language of Croccc & Scoccimarro (2006b, a) or 
r^^' in the notation of Bernardeau et al. (2008). That is, 
P123 can be obtained from the tree-level expression in equa- 
tion (17) by replacing the linear growth D\ implicit in Pq 
as D\ — >• D\{1 + P13/P0). Similarly, -8(14 corresponds to 
a redefinition (or re-summation) of the P2 kernel. In turn, 
B{2z and B222 are leading terms whose corrections appear 
at higher order in the PT series of equation (16). The re- 
summed kernels have well defined properties in terms of tree- 
level quantities and might be the window to an accurate 
description of the non-linear bispectrum at nonlinear scales. 

In Fig. 1 we show the different components in PT 
to the equilateral configurations {upper panels) and to the 



squeezed configurations (lower panels) of the matter bispec- 
trum, respectively B(k,k,k) and B{k,k, Ak) with Ak ~ 
0.012 ft Mpc~^ as a function of k. In the left panels, we 
compare the tree-level contributions Sm with /nl = 100 
{short-dashed, red lines) and B112 {long-dashed, blue lines) 
to the sum of the one-loop corrections P222 + P123 + P123 + 
B114 {dot-dashed, blue lines) present for Gaussian initial con- 
ditions and to the sum of the one-loop corrections due to pri- 
mordial non-Gaussianity, 5(12 + ^^122 + S1I2 + BIjj + -8(^3, 
with fNL = 100 {dotted, red lines). The central and right 
panels compare these sums of one-loop corrections to their 
individual contributions. For the "Gaussian" piece {central 
panels), notice that we plot — Pfls and — S114, implying that 
the overall one-loop correction is the result of a number of 
cancellations similarly to those occurring for the one-loop 
corrections to the matter power spectrum. Analogous con- 
siderations also apply to the "non-Gaussian" one-loop cor- 
rections (right panels) where such cancellations strongly de- 
pend on the triangular configuration. 

To conclude the section, we note that the "order" of 
each correction in the perturbative expansion is defined in 
terms of the power of linear matter density field, j'^' . An 
alternative convention could be given by counting the powers 
of the Gaussian contribution to the curvature perturbations, 
that is <j) in equation (1). Our choice is motivated by the 
standard use in the large-scale structure literature, and by 
the fact that to the n-th perturbative order corresponds the 
well defined rcdshift dependence D^(z). On the other hand, 
we can keep track of the expansion in (j> in terms of the 
nonlinear parameter fNL- 



3 SIMULATIONS 

We utilize a scries of large 1024'' N-body simulations of the 
ACDM cosmology seeded with Gaussian and non-Gaussian 
initial conditions (Desjacques et al. 2009). The box size is 
1600 h~^ Mpc with a force resolution of 0.04 times the moan 
intcr-particlc distance. The (dimensionless) power spectrum 
of the Gaussian part 0(x) of the Bardeen potential is the 
usual power-law A|(A;) = k^ P^{k) / {2n^) = A^{k/ko)"'~^ . 
The non-Gaussianity is of the local form described above. 
We adopt the standard (CMB) convention in which $(x) 
is primordial, and not extrapolated to present epoch. We 
assume h = 0.7, flm = 0.279, fib = 0.0462, Us = 0.96, 
and a normalization of the Gaussian curvature perturba- 
tions = 7.96 X 10~^° at the pivot point feo = 0.02Mpc~\ 
close to the best-fitting values inferred from CMB measure- 
ments (Komatsu et al. 2009). This yields a density fluctu- 
ations amplitude as — 0.81 when the initial conditions are 
Gaussian. Eight sets of three simulations, each of which has 
fNL = 0, ±100, were run with the N-body code gadget 
(Springcl 2005). The same Gaussian random seed field is 
employed in each set of runs so as to minimize the sampling 
variance. The initial particle distribution is generated at red- 
shift Zi = 99 using the Zel'dovich approximation (Zel'Dovich 
1970). 



3.1 Bispectrum estimator and triangle bins 

Let us now introduce the bispectrum estimator B{ki, k2, fcs) 
used in the analysis of the N-body simulations. For a cubic 



Equilateral configurations B{k, k, k) vs. k, f^L = 100, z = 0: 




Squeezed configurations B{Ak,k,k) with Afc = 0.012 /iMpc ^ vs. k, /nl = 100, z = 0: 
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Figure 1. Different components to the PT one-loop prediction for the equilateral configurations B(k, k, k) [upper panel) and the squeezed 
configurations B(Ak, k, k) with Afc ~ O.Ol/iMpc"^ [lower panels) of the matter bispectrum. The left panels show the full prediction 
at one-loop assuming /ml = 100 [black, continuous line), together with the tree-level components Bm [short dashed, red) and B(j^2 
[long-dashed, blue) and the one-loop corrections present for Gaussian initial conditions [dot-dashed, blue) and those depending instead on 
the initial bispectrum and trispectrum [dotted, red). The central panels show the individual terms of the one- loop corrections for Gaussian 
initial conditions, while the right panels show the individual components of the one-loop corrections present only for non-Gaussian initial 
conditions. 



box of volume V, this is given by (Scoccimarro et al. 1998) 

B=^ f d'qi f d'q2 1 d\3 fc(qi23) 5qi ^q,, , (27) 

where Vf = k^ = {2Tr)^ /V is the volume of the fundamen- 
tal cell and where each integration is defined over the bin 
qi £ [ki — Ak/2, ki -\- Ak/2] centered at ki and of size Afc 
equal to a multiple of the fundamental frequency fc/. In our 
case we assume a bin size Ak = 3fc/. The Dirac delta func- 
tion 5o(qi23) ensures that the wavenumbers qi, q2 and qs 
indeed form a closed triangle, as imposed by translational 
invariance. The normalization factor Vsiki, fc2, fcs) given by 

Vs(fci,fc2,fc3) = d^qi d^q2 / d^ga 5D(qi23) 

J ki J k.2 ^3 

~ Stt^ k^k2kz Ak^ (28) 

represents the number of fundamental triangular configura- 
tions (labelled by the triplet qi, q2 and qa) that belong to 
the triangular configuration bin defined by the triangle sizes 
fci, k2 and fcs with uncertainty Ak. In order to better inter- 
pret the simulation results, we provide as well the expression 
for the variance of the bispectrum associated with this esti- 
mator. At leading order, the variance reads as (Scoccimarro 
et al. 1998) 

AB"(fcl,fc2,fc3) = Sfl^P(fcl)P(fc2)P(fe3), (29) 



where the symmetry factor SB{k\,k2,kz) = 6, 2 or 1 for 
equilateral, isosceles or scalene configurations. This expres- 
sion neglects further corrections depending on the matter 
bispectrum, trispectrum and six-point functions that are re- 
sponsible for correlations between different configurations 
(see Sefusatti et al. 2006). 

When comparing the measured bispectrum configura- 
tions to the theoretical predictions in perturbation theory, 
one should be careful to properly account for the effect of 
the finite size of the triangle bins. As explained above, each 
configuration is defined in terms of the sides of the trian- 
gle with ki being the central value and Afc the uncertainty. 
Since we are assuming Afc = Skf, a typically large number 
of "fundamental" triangles fall into each triangle bin. For 
instance, it is easy to see that, in the case of equilateral con- 
figurations, the bin defined by the central value k will include 
equilateral triangles of side q = k — kf or q — k-\- kf just as 
well as nearly-equilateral triangles with different sides still 
belonging to the fc-bin. What is important here is the fact 
that, in the case of equilateral configurations, we will have 
slightly more triangles of size larger than the fundamental 
equilateral triangle with side q = k, than triangles of smaller 
size. This simply follows from the larger number of modes 
at larger q. 

The correct approach consists in computing the raw PT 
prediction B {qi,q2,q3) and average it over the triangle 



bin defined by fci, fe, fcs and Ak, that is 

S*''(fcl,fe2,fc3) = :^ /dV /rf\2 fo(qi23) X 

5^^(91,92,93), (30) 

where S*** is the value to be compared with the measure- 
ments. This is, however, computationally challenging espe- 
cially in the case of the one-loop corrections to the bispec- 
trum, which usually involve three-dimensional integrations. 
An alternative solution, less rigorous yet reasonable given 
the uncertainties of our measurements, consists in defining 
the following effective values ki for the wavenumbers ki char- 
acterizing the triangle, 

ki = ji- I d^qi / d^q2 I d^qs <5D(qi23) qi , (31) 

so that the theoretical prediction which the binned measure- 
ments of the bispectrum must be compared to is 

B*''{ki,k2,k3) = B'"^{h,k2,k3) . (32) 

This procedure improves significantly the agreement be- 
tween theory and simulations, particularly for "squeezed" 
configurations where fci <C ^2 — A;3. 

Here and henceforth, all theoretical predictions will be 
computed in terms of the effective triangle fci, k2 and fes 
as defined above. Furthermore, when the bispectrum is ex- 
pressed as a function of the angle 9 between two of the three 
wavemodes, it is convenient to introduce an effective angle 
9 given by 

cos9{ki,k2;k3) = d^qi d^q2 5d {qi23) x 

Vb Jki Jk2 J k3 

C08 9{qi,q2;q3) , (33) 

where ^(91,92; 93) is the angle between the vectors qi and 
q2. This expression defines the effective angle as a weighted 
average of the angles corresponding to the "fundamental" 
triangles falling in a given bin. These are limited by the tri- 
angle inequalities ^ qi + q2 and qs ^ |gi — g2|. In the 
figures, the (juantitios measured in the N-body simulations 
will be plotted as a function of 9, while the theoretical ex- 
pectations will be the raw PT predictions. 

3.2 The power spectrum 

To facilitate the comparison between different statistics and 
help interpreting the bispectrum measurements of the next 
section, we will first present measurements of the matter 
power spectrum, highlighting the effects due to primordial 
non-Gaussianity and their description in PT. Similar results 
can be found in Desjacques ct al. (2009); Grossi et al. (2008); 
Pillepich et al. (2010); Bartolo et al. (2009). 

In the two upper rows of Fig. 2 we show the matter 
power spectrum measured in simulations of Gaussian ini- 
tial conditions, as well as the linear {dashed lines) and one- 
loop (continuous lines) predictions in PT. In addition, in 
the second row, displaying the ratio of the Gaussian power 
spectrum with a smooth {i.e. no- wiggles) linear power spec- 
trum, we show the nonlinear power spectrum obtained with 
the halof it code of Smith et al. (2003) {thin, green line) 
and the predictions in Rcnormalized Perturbation Theory 
(RPT) (Crocce & Scoccimarro 2006b,a, 2008). The various 



columns correspond, from left to right, to the redshift z = 0, 
1 and 2, respectively. The well-known failure of one-loop PT 
to describe the matter power spectrum at mildly nonlinear 
scales and low redshift is quite apparent (see Crocce & Scoc- 
cimarro 2008; Jeong & Komatsu 2006, for a recent compar- 
ison with simulations and for a comparison at high redshift, 
respectively). The halof it prediction is significantly better 
at low redshift, while it shows a discrepancy of the same or- 
der at z = 2. On the other hand, the RPT prescription pro- 
vides very good predictions (within 1%) up to 0.23ftMpc~^ 
at redshift zero, and over the whole range we consider at 
redshift z = 1 and 2. The slight discrepancy at z = 2, of the 
order of 0.5%, not present at z = 1, might be perhaps be 
explained in terms of transients from the initial conditions 
(Scoccimarro 1998; Crocce et al. 2006), despite the relatively 
high redshift {z = 99) Eissumed for the simulations. 

In the third row, we show the ratio between the matter 
power spectrum extracted from the /jvl = 100 and Gaus- 
sian simulations. The plots for z = and z — 2 repro- 
duce Fig. 3 in Desjacques et al. (2009). Finally, the last row 
shows the difference between the two cases, i.e. P{k; Jnl = 
100) — P{k;fNL = 0). In all these plots, the ratio and the 
difference measured in the simulations are computed for 
each realization and then averaged over the eight realiza- 
tions available. At redshift zero, the one-loop correction P12 
reproduces qualitatively the effect due to primordial non- 
Gaussianity, but it breaks down at relatively large scales, 
k ~ 0.2/iMpc^^, maybe suggesting the need for higher or- 
der corrections. An extension of perturbation theory such 
as the time-renormalization group approach (Matarrese & 
Pietroni 2007; Pietroni 2008) seems to improve only in part 
the agreement between theory and simulations beyond this 
scale (see Fig. 4 in Bartolo et al. 2009). 

In each realization of the initial conditions with /^/^ = 
0, ±100, we also measured the combination [P{k;fNL = 
+100) + P{k; fNL = -100) - 2P{k; fNL = 0)]/2. In the PT 
framework, the result is expected to be the sum of all the 
corrections depending on even powers of /at l ■ At the lowest 
order however, these are given by two-loop contributions 
which we ignore in this work. Nevertheless, we find that in 
the range of scale considered here and for Jnl ~ 100, such 
terms represent an effect of the order of 10^* relative to the 
power spectrum for Gaussian initial conditions. 



4 RESULTS 

We now present the measurements of the matter bispectrum 
with Gaussian and non-Gaussian initial conditions together 
with one-loop PT predictions. In the figures, we will often 
denote these quantities as Bq and Bng, where the "G" 
and "NG" subscripts refer to the initial conditions. In the 
Gaussian case moreover, we will also perform a comparison 
between the measurements and the fitting formula of Scoc- 
cimarro & Couchman (2001). 

To assess the agreement between PT and N-body mea- 
surements as a function of scale and triangle shape, we will 
consider five sets of configurations. We will present results 
as a function of k for equilateral configurations B{k,k,k), 
isosceles configurations B{2k, 2k, k) as well as increasingly 
"squeezed" configurations B{k, k, Ak) with fixed Ak. To fur- 
ther explore the shape dependence, we will also show the 
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Figure 2. Measurements of tlie matter power spectrum, P{k), as a function of k. We sliow, from top to bottom, the power spectrum B 
{first row) and its ratio to the no-wiggle, linear prediction {second row) for Gaussian initial conditions, the ratio P{fNL = +100)/P{/jvL = 
0) {third row) and the difference P{fNL = +100) — P{fNL = 0) {last row). Different columns correspond to redshifts z = 0, 1 and 2. 
Short-dashed, black line indicate the tree-level FT predictions while continuous, black lines the one-loop ones. In addition, on the second 
row we include the RPT prediction of Greece & Scoccimarro (2006b) (at the two-loop approximation, dot-dashed, red line) and the 
prediction from the code halof it of Smith et al. (2003) {long-dashed, green line). 



result of measuring the matter bispectrum for two sets of 
generic configurations for which the magnitude of two sides 
of the triangle (fci and k2) is fixed while the angle 9 between 
them is varied. 

For each of these sets, in each of the following figures, 
the upper two panels show measurements of the matter bis- 
pectrum B, or the reduced bispectrum Q, see equation (34) 
below, for Gaussian initial conditions, as well as the ratio 



to the corresponding tree-level expression in FT where the 
acoustic oscillations are removed by means of the smooth 
transfer function of Eisenstein & Ifu (1998). Recall that 
there is no "linear" matter bispectrum for Gaussian initial 
conditions (but there is an initial bispectrum in the presence 
of primordial non-Gaussianity). For sake of comparison, we 
take the tree-level prediction as a reference since it is most 



directly related to the linear bispectrum, which is generically 

The last three rows in the plots focus on the effect of 
primordial non-Gaussianity. We show, in particular, the ra- 
tio 

B{fNL = WO)/ B{fNL = 0) {third row) , 
the difference 

B{fNL = 100) - B(/ivL = 0) {fouHh row) 

with respect to the Gaussian case, and the combination 

[B{fNL = +100) + B{fNL = -100) 
-2 B{fNL = 0)]/2 {last row) 

to highlight the effects proportional to f^L- In ^-U cases, the 
N-body results indicate the mean over eight realizations of 
the specific combination (ratio, difference, etc.) performed 
with Gaussian and non-Gaussian initial conditions drawn 
from the same random seed field (p (see section 3). In this 
way, we can study the effect of non-Gaussianity without the 
additional sampling variance affecting, for instance, the dif- 
ference Bng — Bg obtained as the difference between the 
mean Bng arid the mean Bg over the eight realizations. Fi- 
nally, the three columns correspond to the results at redshift 
z = 0, 1 and 2. In all the plots, the numerical results are 
compared to the tree-level {short-dashed, black lines) and 
one- loop predictions {continuous, black lines) in PT. 

In Fig. 3, we show the matter bispectrum B{k, k, k) for 
equilateral configurations. As can be seen, non-linearities are 
particularly severe, consisting in a almost ~ 300% correction 
relative to the trcc-lcvol prediction for k ~ 0.2/iMpc~^ and 
z = for instance. The bispectrum measured from a total 
simulation volume of ~ 33 Gpc'^ presents errors of the 
order of 10% at this scale for equilateral configurations. No- 
tice that this specific triangle shape suffers, unlike other con- 
figurations close in shape and scale, from a relatively large 
variance (up to a factor of six). This effect originates partly 
from the symmetry factor ss in equation (29), and from the 
large contribution of higher-order correlation functions to 
the bispectrum variance. 

The one-loop prediction appears to be well within our 
errors up to A; ~ 0.15 h Mpc~^ and describes reasonably well 
the behavior at smaller scales. For k < 0.15/i.Mpc~'^, the 
one-loop prediction behaves better than the fitting formula 
of Scoccimarro & Couchman (2001) (in the plots SCOl), 
which under-predicts the data points at mildly non-linear 
scales. This ~ 20% discrepancy, unsurprising given the size 
of the simulation box used for the fit (240 Mpc), has al- 
ready been noted by Pan et al. (2007). It should be remarked 
that the SCOl formula aimed at describing the nonlinear 
bispectrum at smaller scales, particulalrly for weak lensing 
applications, and did not addressed specifically the issue of 
the acoustic features. Pan et al. (2007) also proposed a phe- 
nomenological model for the matter bispectrum based on a 
rescaling argument similar to the one explored in Hamilton 
et al. (1991); Peacock & Dodds (1996). We also compared 
this prescription to our measurements and find that it agrees 
better than the fitting function of Scoccimarro & Couchman 
(2001). However, the rescaling induces an large and unphys- 
ical shift in the acoustic oscillations that should be properly 
accounted for (in Pan et al. 2007, comparisons arc shown 
with simulations of featureless matter power spectra) . 



The third row of Fig. 3 shows the effect of pri- 
mordial non-Gaussianity in terms of the ratio B{fNL ~ 
100)/ B{fNL = 0). It is interesting to notice that the ad- 
ditional non-linear contributions due to non-Gaussian ini- 
tial conditions correspond, for those set of configurations, to 
a ~ 5% correction regardless of redshift. In fact, the con- 
tribution of the initial bispectrum Bo to this effect is al- 
ready subdominant at fc ~ O.l/iMpc"'^ and z = 0, while 
one-loop corrections themselves fail to account for it at 
slightly smaller scales. This is also apparent in the difference 
B{fNL = 100) — B{fNL = 0) which, in the PT picture, arises 
from the one-loop contributions depending on the initial bis- 
pectrum and trispectrum. At redshift zero, these provide an 
accurate description of _B(/jvl — 100) — B{fNL = 0) up to 
k 0.15 /tMpc"^ 

Finally, in the last row we compare the combination 
[B{fNL = -1-100) + B{fNL = -100) - 2 B{fNL = 0)]/2 
to -B112 which, in the one-loop approximation, is the sole 
term depending on the initial trispectrum and, therefore, on 
Inl- This term appears to underestimate by about 50% (at 
best) the simulation results. One should nonetheless keep in 
mind these contributions represent a 0.1% correction to the 
matter bispectrum. 

In Fig. 4, we show the matter bispectrum for the isosce- 
les configurations, B{2k, 2k, k), as a function of k. The shape 
of the triangle is unchanged while its size is rescaled. In this 
series of plots, the relatively smaller variance (with respect 
to that of the equilateral shape) expected from the discus- 
sion above is quite apparent. The error on the mean is of the 
order of 2-3% for most of the isosceles configurations consid- 
ered. These small errors allow a more accurate comparison 
of the measurements with PT predictions. Note that, while 
each triangle now involves two different scales k and 2k, the 
results are shown as a function of the smaller one (fe) solely. 
For Gaussian initial conditions, the one-loop predictions sys- 
tematically overestimates the data points by more than 10% 
at 2 = 0, but the agreement substantially improves at higher 
redshift. By contrast, the accuracy of the fitting formula 
of SCOl is reasonably good for all the scales and redshifts 
considered. As for the effect of primordial non-Gaussianity, 
considerations similar to those made for equilateral config- 
urations also hold for the isosceles shape. 

In Fig. 5 we compute B{k, k, Ak) on triangles one side 
of which is held fixed to the smallest available fc-bin Afe while 
the other two are equal and varying, k is increased smoothly 
such that this configuration, which represents the coupling 
between the scales k and Ak, asymptote to the "squeezed" 
triangle shape. The errors on this highly correlated set of 
configurations are dominated by the large variance of the 
small-scale mode Afe and are typically slightly larger than 
10%. Still, the one-loop approximation for the Gaussian case 
breaks down already around k = 0.15 ^iMpc^^ at redshift 
zero. The SCOl formula shows instead the same discrep- 
ancy noted above around k ~ 0.1/iMpc~^ while it pro- 
vides a better fit to the data at larger fc. At higher red- 
shift however, perturbation theory fares better than the fit- 
ting formula. The limitation of the one-loop prediction for 
the Gaussian case is also apparent in the corrections due 
to non-Gaussianity. However, the theoretical prediction for 
the ratio Bng/ Bg is in remarkable agreement with the data 
(third row of Fig. 5) . Note that the largc-fe limit in this set 
of configurations does not correspond to the more common 



Equilateral configurations B{k,k,k) vs. k, Gaussian initial conditions (/jvi = 0): 





Figure 3. Measurements of the equilateral configurations of the matter bispectrum, B{k,k,k), as a function of k. We show, from 
top to bottom, the matter bispectrum B {first row) and its ratio to the no-wiggle tree-level prediction {second row) for Gaussian initial 
conditions, the ratio B{ffj]^ = -f 100)/-B(/]vl = 0) {third row), the difference B{fi^L = +100) — B{fi^L = 0) {fourth row) and the 
combination [B{fff]^ = -1-100) -I- B{f^]^ = —100) — 2B(/jvl = 0)]/2 {last row). Different columns correspond to redshifts z = 0, 1 and 2. 
Short-dashed, black line indicate the tree-level PT predictions while continuous, black lines the one-loop ones. In addition, on the second 
row we include the fitting formula of Scoccimarro & Couchman (2001) {long-dashed, green lines). 



Isosceles configurations B{k,2k,2k) vs. k, Gaussian initial conditions {/nl = 0): 




Isosceles configurations B{k,2k,2k) vs. k, non- Gaussian initial conditions {/nl = 100): 
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Figure 4. Same as Fig. 



3, but for isosceles configurations, B{2k, 2k, k) as a function of k. 



Squeezed configurations B{Ak, k, k) vs. k, Gaussian initial conditions (fNL = 0): 





Figure 5. Same as Fig. 3, but for squeezed configurations, B{Ak,k,k), with Ak = 3kf ~ 0.012 /iMpc ^ as a function of k. 



"squeezed" limit, obtained fixing two sides of the triangle 
and reducing the third one, so we do not expect an in- 
crease in the non-Gaussian component, since at larger k we 
are probing smaller scales and the suppression due to the 
transfer function is larger. Nevertheless, the non-Gaussian 
corrections are relatively large for this set of configuration, 
ranging from 10 to 30% and growing with redshift. This tri- 
angle shape is, among those we consider, the most directly 
comparable to the measurements of Nishimichi et al. (2009) 
at redshift z — 0.5 and, in particular to the central panels 
of their Fig. 3. Our errors are consistent with theirs, and 
the agreement between our data points and tree-level PT at 
« = 0.5 is also reasonable. 

In the last two figures, we consider generic scalene tri- 
angles for which the length of two sides fci and k2 is held 
fixed while the angle 9 between them (and, therefore, the 
length of the third side) is varied. Such a set of triangular 
configurations is useful to illustrate the shape dependence 
of the bispcctrurri as it includes collapsed, flattened and al- 
most equilateral triangles depending on the choice of fci and 
fc2. To further isolate the shape dependence of the matter 
bispectrum from its scale dependence, it is convenient to 
introduce the reduced bispectrum Q{ki, ^2,^3) defined as 

^ ^ B{ki,k2,k3) , , 

^ - P(fel)P(/C2) + P(/Cl)P(/C3) + P{k2)P{k3) ■ ^ ' 

In the following two figure, we will indeed show the reduced 
bispectrum in the first rows instead of the bispectrum itself. 
Notice that the one-loop predictions for the reduced bispec- 
trum are computed from a proper expansion of the denomi- 
nator in terms of the one-loop expression for the power spec- 
trum (see Sefusatti 2009, for details). The quantities shown 
in the other rows arc the same as before. A second differ- 
ence with the previous plots is the fact that the data points 
are plotted as a function of the effective angle 9 defined in 
equation (33) (see section 3.1). 

In figure 6, we consider the specific case ki = 
0.094 /i Mpc^^ and ^2 = 1.5 ki. For these configurations, 
9 < 0.6 TT implies that all three sides are larger than 
0.1/iMpc~^. On these scales, the agreement of the one- 
loop predictions with the measurements at z = is poor, 
as is evident from the first plots on the second row. Er- 
rors on the bispectrum mean arc typically of the order of 
3%. At redshift « > 1 however, the theoretical predictions 
fall within the errors. Rather puzzling is, however, the rel- 
atively poor agreement at z = 2, in fact present already 
in the previous plots and perhaps related to small descrep- 
ancy between RPT predictions and simulations in the power 
spectrum case. The prediction for the relative effect of pri- 
mordial non-Gaussianity, which is about 3% at all redshift, 
is in good agreement with the data regardless of the triangle 
shape {third row). The apparent bump shown in these plots 
results from the low values of the "Gaussian" bispectrum 
for nearly equilateral triangles evident from the plots in the 
first row, rather then a non-Gaussian feature. Instead, the 
larger non-Gaussian signal expected for triangles approach- 
ing the squeezed limit is observable in the "difference" plots 
on fourth row for ~ tt^. Notably, the same feature ap- 

^ When the "direction" of the three wavevectors with sum equal 
to zero, i.e. ki H- ki -|- ks = 0, is taken into account it easy to 



pears also in the component -B112 dependent on the initial 
trispectrum Tq {last row). 

Similar results are found for a second set of triangles 
where the two sides are now much closer in size, ki = 
O.UhMpc'^ and fc2 = 0.15/iMpc~^ (sec figure 7). In this 
case however, the configurations are very close to equilat- 
eral for ^ ~ 0.6 TT. As a result, we observe at 2 = the same 
discrepancy between PT and simulations than that seen in 
figure 3 at small scales. This disagreement is also apparent 
in the plot of the reduced bispectrum. The non-Gaussian 
correction typically is of the order of 3%, but it increases 
significantly in the squeezed limit 9 t:. This behavior of 
the linear and nonlinear components due to primordial non- 
Gaussianity is also evident in the fourth row showing the 
difference Bng — Bq- 



5 CONCLUSIONS 

A rather surprising effect of local primordial non- 
Gaussianity on the large scale clustering properties of bi- 
ased objects has been observed in various rmmerical studies 
over the last years (Dalai et al. 2008; Desjacques et al. 2009; 
Pillepich et al. 2010; Grossi et al. 2009). These results at- 
tracted a great deal of attention as they showed that mea- 
surements of the power spectrum of galaxies and quasars 
from current data sets can lead to constraints on the local 
non-Gaussian parameter /jvi comparable to those of CMB 
observations (Slosar et al. 2008; Desjacques & Seljak 2010b). 
Previous work assumed that the main effect of primordial 
non-Gaussianity is limited to an extra contribution to the 
matter and galaxy bispectrum. Still, even under such incor- 
rect but "conservative" assumption, it has been shown that 
future large-volume redshift surveys will reach a sensitiv- 
ity to a non-zero f^L comparable or better than the CMB 
bispectrum (Scoccimarro et al. 2004; Sefusatti & Komatsu 
2007). The inclusion of the non-Gaussian bias in the analysis 
of the galaxy bispectrum or, better, in a combined analysis 
of the power spectrum and bispectrum, is desirable to re- 
liably assess the potentiality of forthcoming surveys of the 
large scale structure. 

As a first step in this direction, we have measured the 
matter bispectrum for the main classes of triangle shape us- 
ing a set of large-volume N-body simulations seeded with 
Gaussian and non-Gaussian initial conditions of the local 
type. We focused on mildly nonlinear scales, 0.02 < A: < 
0.3/iMpc~'^, presented a wide choice of triangular configu- 
rations of different shapes and obtained a determination of 
the bispectrum with an overall errors of the order of 3-4%. 
Of particular interest in this range of scales are the nonlinear 
corrections induced by gravitational instability due to non- 
Gaussian initial conditions as they generate an additional 
non-Gaussian signal on top of the primordial component. 
For a nonlinear parameter Jnl = 100, we found that the 
amplitude of these corrections range from 3-4% for generic 
triangle configurations up to 20-30% for "squeezed" config- 
urations where we expect most of the signal for local non- 
Gaussianity. We quantified these corrections with the aid of 

sec that tlic "squeezed" limit is obtained for — > tt, rather then 

6 — as one might naively think just considering the the triangle 
defined by the wavenumber magnitudes alone. 



B{ki, k2, 0) vs. 6 with fci ~ 0.1 /iMpc ^ and k2 = 1.5fci, Gaussian initial conditions {/nl ~ 0): 




Figure 6. Same as Fig. 3 but for generic configurations B(ki, k2,6) with fci = 0.094 /iMpc ^ and fc2 = 1.5 fei as a function of tlie angle 
6 between fci and fc2. Notice that the first row now shows the reduced bispectrum Q{ki, k2, k^), equation (34), rather than B{ki, k2, kg). 





Figure 7. Same as Fig. 3 but for generic configurations B{ki, k2,d) with ki = 0.14/iMpc ^ and k2 = 0.15/iMpc ^ as a function of 
the angle 8 between fci and fc2- Notice that the first row now shows the reduced bispectrum Q(ki, k2, k-j), equation (34), rather than 
B{k\,k2, fca). 



the ratio and the difference between the non-Gaussian and 
tlie Gaussian bispectrurn. Our set of eight different realiza- 
tions of those models ensure that our results are robust to 
sampling variance. We considered simulations snapshots at 
redshift z — 0, 1 and 2. Overall, we found that the mag- 
nitude of the correction induced by non-Gaussian effects is 
similar regardless the scale and the redshift. This is due to 
a compensation between the primordial component that de- 
creases with time on the one hand, and the contribution 
from nonlinear structure growth that increases with time on 
the other hand. 

We compared our results with the predictions of Eu- 
lerian perturbation theory, both at tree-level and one-loop 
(Sefusatti 2009). As expected, and similarly to what happens 
for Gaussian initial conditions, the trcc-lcvcl approximation 
fails at relatively large scales, k ~ 0.05 - 0.1/iMpc~^, even 
at high redshift. One-loop corrections extend significantly 
the predictive power of perturbation theory down to mildly 
non-linear scales k ^ 0.3/iMpc^^ at redshift z > 1, simi- 
larly to the case of the power spectrum analyzed byjeong 
& Komatsu (2006). They describe, in fact, the matter bis- 
pectrurn measured in simulations at the few percent level, 
with an even better agreement with respect to the "rela- 
tive" effect of primordial non-Gaussianity on the Gaussian 
bispectrurn. Furthermore, they also show a good qualitative 
agreement with simulations at redshift zero. 
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